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Reliable quantum information processing in the face of errors is a major fundamental and techno¬ 
logical challenge. Quantum error correction protects quantum states by encoding a logical quantum 
bit (qubit) in multiple physical qubits. To be compatible with universal fault-tolerant computations, 
it is essential that the states remain encoded at all times and that errors are actively corrected. Here 
we demonstrate such active error correction on a continuously protected qubit using a diamond quan¬ 
tum processor. We encode a logical qubit in three long-lived nuclear spins, repeatedly detect phase 
errors by non-destructive measurements using an ancilla electron spin, and apply corrections on 
the encoded state by real-time feedback. The actively error-corrected qubit is robust against errors 
and multiple rounds of error correction prevent errors from accumulating. Moreover, by correcting 
correlated phase errors naturally induced by the environment, we demonstrate that encoded quan¬ 
tum superposition states are preserved beyond the dephasing time of the best physical qubit used 
in the encoding. These results establish a powerful platform for the fundamental investigation of 
error correction under different types of noise and mark an important step towards fault-tolerant 
quantum information processing. 


Large-scale quantum information processing requires the correction of errors during computations. In quantum 
error correction a logical quantum bit (qubit) is encoded in a subspace of multiple physical qubits so that errors 
can be actively corrected without affecting the encoded information. A promising way to correct errors in encoded 
quantum states is to perform feedback based on multi-qubit measurements known as stabilizer measurements [1-3] 
(see Fig. la). These measurements are performed non-destructively using extra qubits (ancillas) and are frequently 
repeated to detect errors before they accumulate. The measurement outcomes are then processed in classical logic 
that identifies the error syndrome, and, in order to enable universal computations [ 1 ], active feedback is applied to the 
encoded system to correct errors where needed. The key experimental challenge is to perform such complete error- 
correction cycles including non-destructive stabilizer measurements and real-time feedback well within the coherence 
time. 

Quantum-error-correction protocols have been explored across a range of platforms [4-14]. Pioneering experiments 
bypassed stabilizer measurements by reversing the encoding to correct errors, thus leaving the quantum state un¬ 
protected [5-11]. Recent breakthroughs have enabled the use of stabilizer measurements to passively track errors in 
quantum states and retrieve stored information afterwards through post processing [12-15]. Here we realize complete 
rounds of active quantum error correction on a continuously encoded qubit by exploiting newly-developed stabilizer 
measurements based on an electron spin ancilla with high-fidelity non-demolition readout, encoding in long-lived 
nuclear spins and real-time correction of errors through fast classical logic. 

The three-qubit code considered here corrects a single phase error on any one of the physical qubits. The logical qubit 
is encoded as =alO)L+/?ll)L with 10 )l = ( 1 +X)il+X) 2 l+X )3 + 1 -X)il-X) 2 l-X) 3 )/A 11 )l = ( 1 +X)il+X) 2 l+X) 3 - 
l-X)ij-X)2l-X)3)/-\/2 and ] ± X) = (jO) ± \\))/^/2. The logical qubit operators are Xp = X 1 I 2 I 3 , Yp = Y 1 Z 2 Z 3 and 
Zl = Z 1 Z 2 Z 3 (or their permutations). Errors are detected by measuring the stabilizer generators X 1 X 2 I 3 and I 1 X 2 X 3 
via an ancilla. For an uncorrupted state both measurements yield outcome -|-1 (no error), but for a phase error (a Z 
operation) on just one of the qubits the two measurements give a unique syndrome of -1 outcomes that identifies the 
error. 


RESULTS 

Stabilizer measurements and real-time feedback. Our qubits are three nuclear spins (/ = 1/2, 1.1% 
abundance) surrounding a single nitrogen-vacancy (NV) centre in diamond, whose electronic spin we use as ancilla 
{S = 1; lO)a : ms = 0 and [Ijo : ms = —1) (Fig. lb). At 4 K, the ancilla combines fast control [16], optical single-shot 
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readout [17] and long coherence times [18] (> 25 ms, Methods). We use relatively remote nuclear qubits (coupling 
to the ancilla 20-50 kHz) that are robust against optical excitation of the ancilla and design decoherence-protected 
gates to control them [9,19] (Methods). All three qubits show long dephasing times T2 (Fig. Ic) with the dominant 
natural errors being phase errors. 

The key challenge for implementing stabilizer measurements in this system is that the ancilla-qubit interaction 
is always present: imperfect knowledge of the ancilla state during or after readout dephases the qubits [20-22]. To 
minimize this dephasing, we implement quantum non-demolition measurements of the ancilla by resonant optical 
excitation of |0)a and stopping the excitation within 2 fis upon photon detection (outcome |0)a) to minimize spin 
flips in the optically excited state [23] (Methods). The resulting readout fidelities are Fq = 0.890(4) for |0)o and 
Fi = 0.988(2) for |l)a (average: F = 0.939(2)). Crucially, the post-measurement fidelity after correctly assigning |0)a 
is 0.992, demonstrating the desired non-demolition character. 

To benchmark the stabilizer measurements and real-time feedback, we deterministically entangle two qubits by 
projecting into a Bell state, i.e. a simultaneous eigenstate of XX and ZZ [21,24,25]. First, the qubits are initialized 
in |00), an eigenstate of ZZ, with fidelity 0.910(6). Then a XX stabilizer measurement projects the qubits onto one 
of two Bell states (Fig. Id). We interpret the -1 outcome as an error in the desired state and correct it through 
feedback before performing two-qubit tomography. The deterministically generated entangled state, with fidelity 
F = 0.824(7) (Fig. le), demonstrates the non-destructive nature of the measurement; coherence within the subspaces 
is maintained throughout the measurement and feedback cycle. The complete cycle can be repeated up to 6 times 
within the shortest qubit T^. 

Active quantum error correction on a logical qubit. We now turn to quantum error correction by stabilizer 
measurements. The logical qubit is encoded by taking an arbitrary state \'ip)a = a|0)a -|-/3|l)a prepared on the ancilla 
to the three-qubit state |V')l = Q;|0)l + I3\1)l (Fig. 2a). We characterize the encoding by preparing six basis states 
|0)l, |1)l, I ± X)l = (|0)z, ± |1)l)/-\/2 and | ± Y)l = (|0)l ± i|1)l)/-\/2 and performing three-qubit state tomography. 
The fidelities with the ideal states confirm successful encoding and genuine three-qubit entanglement (Fig. 2b). 

We first investigate the recovery of arbitrary logical qubit states from phase errors. To emulate a general process 
causing dephasing, uncorrelated incoherent errors are applied with variable probability to each physical qubit 
simultaneously (Fig. 3a); for each qubit the error process is E{p) = (1 —pe)lpl + PeZpZ, with p the single-qubit 
density matrix. By controllably applying such errors we characterize the effectiveness of the error correction for any 
process causing uncorrelated errors with equal probability to the qubits. We then measure the stabilizers X1X2I3, and 
I1X2X3, identify potential errors and correct them through feedback. The probabilities to obtain the four different 
error syndromes (inset in Fig. 3b) show the expected symmetry around Pe = 0.5 and match the theoretical prediction 
based on the errors present in the initial states (Fig. 2b) and the average ancilla readout fidelity. 

The protection of the logical qubit is characterized by the process fidelity with the identity (Fig. 3b) (Methods). 
We quantitatively analyse the results by fitting to wFqec + (1 ~ w)Fiinear, where FQEciPe) and FuneariPe) are the 
theoretical curves with and without error correction (rc = 1 indicates ideal robustness against applied single-qubit 
errors). When no error correction is applied we observe the expected linear dependence on the error probability: 
w « 0. In contrast, with quantum error correction w is 0.81(3), and a non-linear curve shape that is characteristic 
for robustness against single-qubit errors is obtained. This result demonstrates that the entropy associated to the 
applied errors is successfully removed from the system. 

Comparisons to an unencoded qubit and the logical qubit without error correction reveal that adding quantum 
error correction on top of a computation does not yet provide a net improvement (Fig. 3b), due to additional 
errors introduced by the initialization, encoding and stabilizer measurements (total of 13 two-qubit gates, 488 ancilla 
refocusing pulses and 6 ancilla readouts/resets). To isolate the errors due to the stabilizer measurements, we compare 
the error-corrected logical qubit to the logical qubit left idle. We further optimize the error correction, by assigning 
the ancilla state with the best readout fidelity (|l)a, Fi = 0.988(2)) to the most likely error syndrome (-1-1, -1-1 - no 
error, inset Fig. 3b), instead of averaging over all assignments as in Fig. 3b. With this improvement, error correction 
outperforms idling for a range of Pe (Fig. 3c); once the logical qubit is encoded, quantum error correction can be 
beneficial. 

Multiple rounds of active error correction. Because a complete round of error correction (2.99 ms) fits 
well within the dephasing time of the physical qubits, we can concatenate multiple rounds to improve the coherence 
of continuously encoded quantum superpositions by preventing the accumulation of errors (Fig. 4a). Three new 
elements are introduced. First, the total error probability Pe is distributed over n rounds, so that the error probability 
per round is p„ = (1 — >/! — 2pe)/2 (Methods). This error model corresponds to errors occurring incoherently, for 
example with a constant rate in time. Second, to investigate dephasing we consider only the protection of the two states 
I ± X)l = I ± X, ±X, ±X) (i.e. a classical bit stored in the phase of a quantum superposition). Third, we exploit the 
intrinsic robustness of the logical qubit to single Z errors by redefining Xl = (A1/2/3-I-/1A2/3-I-/1/2A3 — AiA2A3)/2 
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which is equivalent to performing a round of error correction by majority voting at the end of the experiment [13,14]. 

For a single round of error correction (majority vote only) the average fidelity is higher than for an unencoded 
qubit for any pe (Fig. 4b); adding more (identical) qubits is always beneficial in the repetition code. For pe = 0, 
additional rounds of quantum error correction can only introduce errors, reducing the fidelity (Fig. 4b). For larger p^, 
however, multiple rounds prevent errors from accumulating by dividing the error process in parts that are more likely 
to contain only single errors, which are corrected. In addition, unlike error detection with post-processing [13,14], 
active correction between rounds keeps the probability to obtain 4-1 (no error) high (inset Fig. 4b) and thus maintains 
the advantage of assigning the highest-fidelity ancilla readout to that outcome. As a result, for pe > 0.3, multiple 
rounds outperform a single round of error correction. 

Correcting natural dephasing. Finally, as an example of suppressing errors naturally present in the environment, 
we let the qubits evolve naturally instead of applying errors (Fig. 4c). The error probabilities are now different for 
each qubit because their intrinsic T 2 differ due to their local environments (Fig. Ic). In addition, the errors arise 
from quasistatic detunings due to the fluctuating spin bath, so that they are correlated in time and can evolve 
coherently. Like most environmental errors, such errors might also be suppressed by other methods than quantum 
error correction, for example by polarizing the spin environment [26,27], by refocusing pulses [28] or by isotopic 
purihcation [28-31]. 

The fidelity for the logical qubit with majority voting again starts above the best unencoded qubit, but drops below 
it for larger evolution times (Fig. 4d). Because the error probabilities vary between qubits, an error detected on the 
best qubit becomes more likely to actually correspond to errors on both other qubits and the wrong correction is 
made. An additional round of quantum error correction in the middle of the evolution time now not only prevents 
errors from accumulating by intermediately correcting them, but also interrupts any coherent build-up by projecting 
the errors, thus suppressing them (Fig. 4d). Due to this combination, the logical qubit with error correction yields 
the highest average fidelity for total evolution times between 5 and 19 ms (Fig. 4d). This result demonstrates an 
actively error-corrected logical qubit with an improved dephasing time over the best qubit used in the encoding. 
DISCUSSION 

The presented non-destructive measurements and real-time feedback on encoded quantum states are the key prim¬ 
itives for universal computations on logical qubits and for error-correcting codes that correct both phase and bit-flip 
errors. To reach scalability thresholds, readout and gate fidelities can be further improved, for example through 
optical cavities [32], the implantation [33] or selective growth of defects and isotopes [28,29], and optimal control 
[33]. In a wider perspective, our results can be combined with recently demonstrated entanglement between distant 
NV centres [34,35] to form quantum networks with error-corrected nodes for entanglement purification, quantum 
communication and networked quantum computation [36]. Therefore, these results establish a promising platform to 
experimentally investigate protocols for fault-tolerant quantum information processing under different types of noise 
and error correlations in diverse settings. 


METHODS 

Sample and Setup. We use a naturally occurring Nitrogen-Vacancy (NV) in high-purity type Ila chemical-vapour- 
deposition (CVD) grown diamond with a 1.1% natural abundance of and a < 111 > crystal orientation (Element 
Six). To enhance the collection efficiency a solid-immersion lens was fabricated on top of the NV centre [17,37] 
(Fig. lb) and a single-layer aluminium-oxide anti-reflection coating was deposited [34,38]. The sample temperature 
is T ss 4.2 K and a magnetic field of 403.553(3) G is applied along the NV symmetry axis. 

The ancilla NV electron spin is characterized by a Rabi frequency of 4.3 MHz, a dephasing time T 2 =4.6(2) ps, 
a Hahn echo time T 2 = 1.03(3) ms and a longitudinal relaxation time of 0.43(6) s (due to microwave noise and 
laser background). The coherence time of the ancilla under dynamical decoupling exceeds 25 ms and does not limit 
the experiments (Supplementary Fig. 1). We initialize and readout the ancilla through resonant excitation of the 
zero-phonon transitions of the NV centre (Supplementary Fig. 2). Prior to every experiment the ^"^N nuclear spin is 
initialized by measurement with a fidelity of Fjv = 0.94(3) in mi = —1 [17]. No external electric fields are applied: 
the gates in Fig. lb are grounded. 

Nuclear spin qubit control. The hyperfine interactions for the three nuclear spins are estimated by dynamical 
decoupling spectroscopy [9] (Supplementary Table 1). Building on previous gate designs [9], nuclear gates are realized 
by applying sequences of 7r-pulses on the electron spin of the form (t — it — 2t — it — The number of pulses 

N sets the rotation angle. The inter-pulse delay 2 t determines which qubit is controlled and whether the rotation is 
conditional on the ancilla state. In contrast to previous work [9] we allow the gates to be detuned, providing greater 
flexibility to optimize r and N for gate selectivity and minimal discretization errors. The gate parameters are listed 
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in Supplementary Table 1 and 2. 

The nuclear spins are initialized by swapping with the ancilla electron spin (Supplementary Fig. 3) and are read 
out by mapping the required correlation to the ancilla before reading it out (Supplementary Fig. 4). To obtain 
best estimates for the actual states, the results are corrected for the hdelity of the gates used in the final readout 
(tomography) (details in Supplementary Note 3). Uncorrected data is shown in Supplementary Fig. 11. 

Feedback. Real-time feedback is implemented through a programmable microprocessor (ADwin Pro II) that 
controls the experimental sequence (Supplementary Fig. 5). We exploit feedback in four different ways. First, detected 
phase errors are corrected directly after the stabilizer measurements. Note that analysing errors over multiple rounds 
[14] would additionally enable real-time correction of ancilla readout errors, but that this is not implemented here. 
Second, dependent on the ancilla measurement outcome the qubits pick up a deterministic phase shift due to the 
hyperhne interaction, which is corrected in the same way. Third, for an odd number of -1-1 outcomes the operations 
in the stabilizer measurements imprint a bit flip on the logical qubit, which we correct by transforming the logical 
qubit basis in real time. Fourth, to start each measurement sequence with the ancilla in |0)o it is flipped back to |0)a 
when the previous measurement returned |l)a- 

Importantly, we perform real-time feedback either by adapting the qubit bases for all subsequent gates and mea¬ 
surements (for correcting Z errors and for the logical qubit) or by absorbing the feedback operations into the next 
gate acting on the same qubit (for the ancilla). Therefore the physical control sequence is directly adapted based on 
the measurement outcomes without introducing unnecessary gate operations that would cause additional errors. In 
the circuit diagrams we sometimes display the gates for the feedback separately for clarity. 

Quantum error correction analysis. The process hdelity with the identity is given by Fp = (Fq -I- Fi -|- F+x + 
F-x + F+y + F-y — 2)/4, with Fa = {a\pa\oi) the six hdelities of the hnal states pa with the ideal states |a)z,. The 
results of Fig. 3 are analysed by htting to wFqEciPe) + (1 - w)Fiinear{Pe), with FqEciPe) = O + A{1 - ipl + 2pl) 
and FiineariPe) = O + A{1 — Pe). A and O account for the experimental hdelities (Supplementary Note I). 

The state hdelities for multiple rounds of error correction and incoherent errors (Fig. 4b) are htted to the same 
equation using FgEciPe) = y (1 + ^(1 ~ -|- 4p^)"), with n the number of rounds, Pn = ^ ~ 2pe the error per 

round, and FuneariPe) = |(1 + ^(1 ~ 2pe))- The error per round is obtained as follows. An error process with 
total error probability (pe) reduces the expectation value by a factor of (I — 2pe). For incoherent errors, a process 
can be divided in n equal rounds using (1 — 2pe) = (1 — 2p„)”, which results in = (I — — 2pe)j2 (for pe < 0.5). 

In Fig. 3c and Fig. 4b, A depends on the error-probability Pe, because we optimize the effective readout hdelity by 
associating the most likely error syndrome to the best ancilla readout (Supplementary Note 1). See Supplementary 
Information for further details on all theoretical analysis, including the error syndrome probabilities and numerical 
simulations of Fig. 4d. 
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Figure 1 : Quantum error correction and implementation of stabilizer measurements, a, A quantum state is 
encoded in a logical qubit consisting of three physical qubits. Errors inevitably occur, for example during computations. An 
ancilla is used to repeatedly perform measurements that detect errors. Errors are corrected through classical logic and 
feedback, while the quantum state remains coherent and encoded, b, Device: CVD-grown single-crystal diamond with a 
solid-immersion lens [37] and on-chip lines for microwave control. Ancilla: the optically addressable electronic spin of a 
nitrogen vacancy (NV) centre. Qubits: three nuclear spins that are controlled and measured through the hyperfine 
coupling to the ancilla (Methods), c, Free induction decay (Ramsey) experiments. Gaussian fits yield dephasing times T 2 = 
= 12.0(9), 9.1(6) and 18.2(9) ms for qubits 1, 2 and 3, respectively, d, Deterministic entanglement of two qubits by XX 
stabilizer measurement and feedback. The ±x gates are n/2 rotations around x with the sign controlled by the ancilla state. 
The final operations reset the ancilla and account for an additional flip for the -Fl outcome (Methods), e, State tomography 
of the generated entangled state for qubits 2 and 3. The fidelity with the ideal state is F = 0.824(7) (see Supplementary Fig. 
6 for other qubit combinations and post-selected results). 
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Figure 2 : Encoding of the logical qubit, a, Encoding an arbitrary quantum state \'(p) = a|0) + /3|1) prepared on the 
ancilla into = a|0)L + /3|1 )l- Successful encoding is heralded by outcome |0)a- b, Characterization of the logical states 
I + X)l, \ +Y}l and |0 )l- Only the logical qubit operators and stabilizers are shown (see Supplementary Fig. 7 for complete 
tomography of all 6 logical basis states). The fidelities with the ideal three-qubit states are F = 0.810(5), 0.759(5) and 
0.739(5), respectively, demonstrating three-qubit entanglement [10]. The logical state fidelities are 

F+x = (1 -E {Xl))/2 — 0.916(6), F+y = (1 -E {Yl))/2 = 0.822(7) and Fo = (1 -E (Z_l))/ 2 = 0.813(9). Ideally, all the encoded 
states are -El eigenstates of the stabilizers X1X2I3 and I1X2X3. The fidelity to this code space, 

Fs = (1 + (X 1 X 2 / 3 ) -E {I1X2X3) + (Xi/ 2 X 3 ))/ 4 , is 0.839(3) averaged over all states and gives the probability that the 
starting state is free of detectable errors. 
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Figure 3 : Active quantum error correction by stabilizer measurements, a, All qubits are simultaneously subjected 
to uncorrelated phase errors E with probability pe- Errors are detected by measuring X\X 2 l 3 and I 1 X 2 X 3 and subsequently 
corrected by Z operations through feedback. Finally we measure the process fidelity with the identity, b, Process fidelities 
for: an unencoded qubit, the logical qubit without stabilizer measurements, the error-corrected logical qubit, and the logical 
qubit without feedback (i.e. errors are detected but not corrected). We average over the logical qubit permutations, e.g. 

Xl = X 1 I 2 I 3 , I 1 X 2 I 3 and I 1 I 2 X 3 , and the four ways to assign the ancilla states to the error syndromes (see Supplementary 
Fig. 8 for individual curves). Inset: probabilities for the error syndromes with theoretically predicted curves based on the 
state tomography in Fig. 2b (Supplementary Note 2). c, Comparison between the error-corrected logical qubit and the logical 
qubit with the stabilizer measurements replaced by an equivalent idle time (2.99 ms). Compared to b, the effective readout 
fidelity is optimized by associating syndrome -1-1,-|-1 (no error) to obtaining |l)a for both stabilizer measurements. Curves in 
b and c are fits described in the Methods. 
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Figure 4 : Extending coherence by active quantum error correction, a, Three rounds of error correction on a 
logical qubit. The first two rounds of quantum error correction use stabilizer measurements and feedback. The hnal round is 
implemented by majority voting, b, Average logical state fidelity for | + X)l and | — X)l as a function of total error 
probability pe for n = 1, 2 and 3 rounds of error correction compared to an unencoded qubit. The errors per round En occur 
with probability p„. Inset: probabilities that no error is detected (n = 3). The similarity of the results for rounds A and B 
confirms that errors are corrected in between rounds, c, Correcting natural dephasing. The storage time is defined from the 
end of the encoding until the start of the final measurements, d, Dephasing of the logical qubit: without stabilizer 
measurements, with quantum error correction and without feedback, compared to the best unencoded qubit. The dashed lines 
indicate the times between which the actively error-corrected logical qubit gives the highest fidelity. The data without 
feedback (detecting errors without correcting) isolates the suppression of coherently evolving errors by projecting them. For 
long times, applying error correction lowers the fidelity because the stabilizer measurements extract no useful information 
about errors, but nevertheless preferentially suppress evolutions that result in phase errors at the end of the sequence (See 
Supplementary Fig. 10 for a detailed analysis). See Supplementary Fig. 9 for error syndrome probabilities. Curves in b and d 
are fits described in the Methods and Supplementary Notes 1&2. 
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SUPPLEMENTARY INFORMATION 

Supplementary Note 1 : Theoretical analysis: state and process fidelities 

For ideal error correction the process fidelity to the identity as a function of error probability Pe for a single round 
of quantum error correction (QEC) is 

FQEc{Pe)=0 + A{l-ipl + 2pl). (1) 

The offset O and amplitude A account for the finite experimental state fidelities. Note that the value at Pe = 0-5j 
FqEciPe = 0-5) = O + A/2, is determined by the fidelity of the logical states |0 )l and |1 )l, which are insensitive to 
phase errors. Without error correction a linear function 

Flinear = O + A{1 — Pe) (2) 

is expected. The experimental data can be fitted to a weighted sum of the two functions by: 

Fp = wFqeC + (1 — w)Fiinear- (3) 

The shape of the curve is set by w, which gives the relative weights of the equations for ideal error correction and for 
no error correction. 


Assignment of ancilla states to the error syndrome: effective measurement fidelity Fm 

In our experiment, the ancilla readout fidelity is asymmetric: |l)o has a higher readout fidelity {Fi = 0.988(2)) 
than |0)a (Fb = 0.890(4)). The effective measurement fidelity for error correction Fm is therefore determined by the 
probabilities to obtain |0)a or |l)a, which depend on the assignment of the ancilla states (|0)a or |1 )q) to each stabilizer 
measurement outcome (+1 or -1) and the probabilities for different errors to occur. There are four different ways to 
assign the ancilla states to the error syndromes: the +1,+1 outcome (no error) can be set to result in {|0 )q, |0)a}, 
{|0)a, |l)a}, {|l)a, |0)a} Or {|l)a, |l)a}- The probability to obtain outcome +1,+1 (no error) ideally is 1 — 3pe + 3pg, 
while the probability to detect an error on a given qubit is Pe — Pe- With these probabilities, we obtain the effective 
QEC measurement fidelity as function of error probability: 

Fm = F(o)( 1 - 3pe + Spl) + (F(i) + F^^) + F^^^){pe - pi) = F^°^ + (F^^) + F^^) + F^^) - 5F^°^){pe - pi), (4) 

with F(°\ F^^), F^^) and F^^\ the readout fidelities for 0 errors, an error on qubit 1, an error on qubit 2 and an 
error on qubit 3, respectively. Eor example, for assignment {|l)a, |l)a} to stabilizer outcomes +I,+1 (no error), these 
readout fidelities are 


F(o) = 

F(i) = f(3) = FiFo, 

F(2) = F^, 

In a similar way, the fidelities for the other three assignments can be calculated. 

Finally, if we assume that an erroneous ancilla readout decoheres the logical state, the dependence of the 
effective readout fidelity on pe can be taken into account by setting: 

A = A'Fm (5) 

in Eqs. dMH for the process fidelity, with A' a constant. 


Fitting of Figs. 3b, 3c and[^ 

In Fig. 3b, the ancilla readout is symmetrized by averaging over all four assignments, so that Fm equals the 
average readout fidelity 0.939(2) and is independent of Pe- We can therefore simply fit the data in Fig. 3b to Eq. 
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with A constant. We find w = 0.81(3), corresponding to an average probability to successfully correct single-qubit 
errors of (P„) = |(w -I- 2) = 0.94(1) [1]. We obtain A = 0.557(2) and O = 0.086(1). For the unencoded qubit, the 
encoded qubit without stabilizer measurements, and the encoded qubit without feedback, we find a linear function 
and w « 0 ((P„} « 2/3) as expected without error correction (exact values: w = —0.06(3),—0.03(3) and —0.07(3), 
A = 0.882(4), 0.734(3) and 0.543(2), and O = 0.019(3), 0.051(3) and 0.092(1), respectively for the three cases). In 
Fig.0 the separate process fidelities for the different assignments are shown. Switching between assignments is done 
by adding or omitting a 7r-pulse before the ancilla readout. 

In Fig. 3c, we assign the ancilla state |l)a to the -1-1 outcome for all stabilizer measurements. This assign¬ 
ment is optimal because it associates the best readout fidelity with the most likely outcome: -1-1,-1-1 (no error, inset 
in Fig. 3b). We fit the data in Fig. 3c and Fig. [^to Eq. § with A now error-dependent according to Eq. and obtain 
w = 0.8(1), corresponding to (P„) = 0.93(3) (A' = 0.666(8) and O = 0.038(6)). The values for w and (P„) are in 
good agreement with the result of Fig. 3b, indicating that the treatment in Eqs. |4]^^is accurate. 


Multiple rounds of error correction (incoherent errors), Fig. fb 

For multiple rounds of QEC with incoherent errors and with the total error with probability Pe equally distributed 
over n rounds, the error-probability per round is Pn = ^ (1 ~ ^1 — 2pe), for pe < 0.5. Ideally, the (average) state 
fidelity is then described by: 

F=^[l + il-6pl+Aplr]. (6) 

As before we fit the data to a weighted sum of the equations for ideal error correction and for a linear error-dependence 
(no error correction). We use the optimal ancilla state assignment {FuiPe) from Eq. |^. For two rounds of error 
correction we obtain 

F2 = izc[l + A'Fm{ 1 - Gpl + ^plf] + i(l - u;)[l + A'Fm(1 - 2p,)], (7) 

giving w = 0.66(4) and A' = 0.850(9). For three rounds it becomes 

P 3 = ~ 6P3 + 4p|)^] -I- -(1 — w))! -I- A'F^(1 — 2pe)], (8) 

giving w = 0.71(2) and A' = 0.810(5). Importantly, the data for multiple rounds cannot be accurately described by 
the expected shape for a single round of error correction (Eq. [^. 


Naturally occurring decoherence (coherent errors), Fig. fd 

The experiments for the best unencoded qubit and logical qubit (majority vote only) are fitted to a general expo¬ 
nentially decaying function: 

F= ^(l + Ae-(‘/^F). (9) 

Here, we obtain for the best qubit: T = 17.3(2) ms and n = 2.09(7), while for the encoded qubit with majority voting 
we obtain: T = 13.7(1) ms and n = 2.37(8). 

To understand the interplay of quantum error correction and the projection of errors in the experiments with 
stabilizer measurements at half the free evolution time in Fig. 4d we turn to numerical Monte Carlo simulations, see 
Fig. [ 1 ^ for details and results. 


Supplementary Note 2 : Theoretical analysis: error probabilities 


The probability to detect no error ) is the sum of the probability to have no error (no qubits flipped) or three 
errors (all qubits flipped) and is described by: 

p(0) = (1 _pW)(l _p(3)^ + p[l}^p[^^p[^^, 


( 10 ) 
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where is the error probability for qubit i. The probability to detect an error on one of the three qubits is 

the probability to have an error on qubit i, or an error on both of the other qubits, which for example for qubit 1 is 
described by: 


=Ptokl -Ptob(l -Ptol) + (11) 

With finite input error probability for qubit i (errors already present in the initially prepared state), the total 
error as function of the applied error probability pe, becomes: 

Ptl = P|n +Pe- 2p[nVe (12) 


Finally we can take imperfect ancilla readout into account and obtain the probability to detect one of the error 
outcomes Pq ^ (i =0 for no detected error) as function of the applied error p^. 


pW ^ p(0)p2 p ^p(l) _ p) + p(2)('l _ p)2 (-^3^ 

= p(l)p^2 ^p(O) p(2))p(l _ p) + p(3)(l _ p)2 (^4) 

pW ^ p{2)p2 p(3))p(i _ p) p p(0)('l _ p)2 (-^5^ 

p0) ^ p(3)p2 p ^p(O) p p(2))p('l _ p) + p(l)(l _ p)2 (10) 


The XX stabilizers in the encoded state tomography (Fig. 2) detect errors present in the encoded state, we 
obtain: 


p(0) =(l(l + Xi,X2,/3)(l + Xi,J2,X3)) 
PW =(^(l-Xi,X2,/3)(l-Xi,/2,X3)) 
p(2) =(l(l-Xi,X2,/3)(l + Xi,/2,X3)) 
P(3) =(^(l + Xi,X2,/3)(l-Xi,/2,X3)) 


i(l + (Xi,X2,/3) + (/l,X2,X3) + (Xi,/2,X3)) = 0.785(2) 
i(l- (Xi,X 2 ,/ 3 ) + (Jl,X 2 ,X 3 )- (Xi,/ 2 ,X 3 )) = 0.060(2) 
i(l- (Xi,X 2 ,/ 3 )- (/i,^ 2 ,^ 3 ) +^ 3 )) = 0.083(2) 
^(1 + {X,,X2,h) - (/l,X2,X3) - (Xi,/2,X3)) = 0.071(2) 


which are uncorrected for qubit readout. These results can be translated to the input errors, as these outcomes refer 
to Eqs. 10?^11 with no additional applied error pe, giving p[^^ = 0.064(2), p[^^ = 0.091(2), = 0.077(2). 


Using these values we estimate the expected total error detection probabilities Pj^\ P^\ Pj^^ and as 

function of applied error probability p^ according to Eqs. HEl The expected error-dependent QEC measurement 
outcomes are shown by the solid lines in the inset of Fig. 3b. 


Error syndrome assignment 

For the different error assignments, the asymmetry in the ancilla readout complicates the error detection curves: 
the QEC measurement fidelity is dependent on the error probability. If, for instance, both stabilizer measurements 


giving -1-1 are assigned to {|l)a, |l)a}, Eqs. mill become: 

p{0) ^ p(0)p2 p (p(l) p p(3))p^(4 _ pg) p p(2)(i _ p^)2 (17) 

pW ^ p(i)p^pQ P p(0)p^(i _ p^) P p(2)Pq(i _ Pq) p p(3)(i _ p^)(i _ Pg) (18) 

p{2) ^ p(2)p2 (p(l) p p(3))pg(l _ p^) p p(0)(l _ p^)2 (ig) 

pi3) ^ p(3)p^p^ p p(0)p^(i _ p^) p p(2)Pq( 1 _ pg) p p(2)(i _ p^)(i _ p^) (20) 


All error detection curves for the four error assignments using similar equations are plotted in Fig. 
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Multiple rounds of error correction, Fig. 4 b 

For multiple rounds we now calculate the average input error from the detection probability for no additional 


applied error (pe = 0). We simplify Eq. 10 


to 




( 21 ) 


and use Eq. [^to obtain the following average input error for round 1: = 0.092(1) and for round 2: 

0.086(1). The resulting curves according to Eq. 17 are shown in the inset of Fig. 4b. 


Supplementary Note 3 : Qubit readout calibration 

To obtain best estimates for the actual states, the results are corrected for the fidelity of the gates used in the final 
readout (tomography). We distinguish between reading out single- two- and three-qubit expectation values. 

For a single qubit i that is initialized and readout immediately, the measured expectation value {Zi) is set by 
the initialization hdelity of the nitrogen spin (F)v = 0.94(3)) and by factors due to the initialization (Cinit.Qi) and 
readout (Cq.) of the qubit. Because the initialization and readout consist of the same set of gates, we assume that 
Cinit,Qi = Cq. for this experiment and obtain: 

{Z,) = F^C^q^, ( 22 ) 

from which a readout correction factor IjCg. can be determined. 

To calibrate the multi-qubit readouts we initialize the three qubits in separable states. For example, for state 
|000), the measured three-qubit expectation value {Z 1 Z 2 Z 3 ) is set by the nitrogen initialization F^, by factors Cinit,Qi 
due to the individual initialization hdelities of the three-qubits and by a factor Cqj_q 2 ^q 3 due to the three-qubit 
readout: 


{Z 1 Z 2 Z 2 ,) — FivCinit,QlCinit,Q2Cinit,Q3CQj,Q2,Q3 t Cq,^^Q2,Q^ 


{Z1Z2Z3) 

Fn Cinit, Q 1 Cinit ,Q2 Cinit, Q 3 


(23) 


This equation assumes that the initialization errors, other than those due to the nitrogen initialization, are uncorre¬ 
lated. The initialization fidelities are obtained using the single-qubit expectation values and single qubit Cq. for the 
corresponding qubit, i.e. for qubit 1 : 


{^ 1 ^ 2 ^ 3 ') — -C/vCjnit.QlCQ^ i Cjnit.Ql — 


{Zihh) 

FnCq^ 


(24) 


with Cqj from Eq. 22 In a similar way, the two-qubit readout is calibrated using two-qubit expectation values of 
two- and three-qubit states. We obtain the following values: 


Cq, = 0.95(1) 
Cq 2 = 0.94(1) 
Cq 3 = 0.95(1) 


Cq,q 2 = 0.94(2) 
Cq,q 3 = 0.88(4) 
CQ 2 Q 3 = 0.90(2) 


Cq.q.q, = 0.92(5) 


which are used to calibrate the final readouts for tomography. Note that the uncertainty in the readout calibration 
potentially creates a small systematic error (a rescaling of all y-axes). Eor this reason we also provide all raw 


(uncalibrated) data for the error correction in Fig. 11 
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SUPPLEMENTARY TABLES 



Qubit 1 

Qubit 2 

Qubit 3 

All (kHz) 

27r-20.6 

27r-36.4 

277-24.4 

Au (kHz) 

27r-43 

27r-25 

277-26 

cJo (kHz) 

27r-431.874(3) 

271-431.994(3) 

277-431.934(3) 

cJi (kHz) 

27r-413.430(3) 

277-469.025(3) 

277-408.303(3) 

r (^s) 

13.616 

4.996 

11.312 

N 

32 

34 

48 

gate time (/rs) 

980 

400 

1086 

T2*, rUs = 0 (ms) 

12.0(9) 

9.1(6) 

18.2(9) 

T 2 , ms = -1 (ms) 

12.8(6) 

9.8(4) 

21(1) 

Ti, rUs = 0 (ms) 

110(10) 

100(10) 

330(30) 


Supplementary Table 1 : Qubit and gate parameters. Ay and A± are the estimated hyperfine interaction components 
parallel and perpendicular to the applied magnetic field, loq and uji are the nuclear precession frequencies for nis = 0 (|0)a) 
and rus = —1 (|l)a). ris half the inter pulse delay, N the number of pulses and gate time the total duration for the 
conditional ±a:-gates. These values vary slightly over the experiment as they are calibrated every ~ 36 hours. T 2 is the 
(natural) dephasing time and Ti the longitudinal relaxation time. 



Fig. 3b (logical qubit with QEC) 

Fig. 4b (three rounds) 

Two-qubit gates 

19 

20 

Ancilla refocussing pulses 

698 

808 

Ancilla read-out and reset 

7 

9 


Supplementary Table 2 : Experimental complexity. Number of operations in the sequence starting from the 
initialization of the qubits. All qubit (^^C) gates are composed of ancilla (NV electron spin) refocussing pulses and the ancilla 
is read-out and reset multiple times. We give values for two examples; a single round of QEC with measurement of 
(Fig. 3b) and three rounds of QEC with measurement of (X1A2A3) (Fig. 4b). 
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SUPPLEMENTARY FIGURES 



Supplementary Figure 1: Coherence of the NV electron spin by dynamical decoupling with 256 pulses. The 

spin is decoupled from the nuclear spin bath by applying a sequence of 256 7r-pulses with alternating phases. The time 
between the pulses is chosen to be a multiple of the Larmor period of the spins. This result shows no significant decay on 
the relevant timescale of our experiments. 
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Supplementary Figure 2: Optical readout and initialization of the ancilla NV electron spin. The electron spin is 
initialized and read out by spin-selective resonant excitation [2]. To initialize or reset the electron spin state we apply a laser 
pulse that excites only the rria = ±1 -f-j- transitions (Reset). Due to spin mixing in the excited state this prepares the 

electron spin in the nia = 0 state (fidelity > 0.98) [2]. To measure the spin state we apply a laser pulse resonant with the 
nia = 0 -H- Ex transition (Readout). Ideally, this results in the detection of 1 or more photons for the rris = 0 state, and no 
detected photons for nia = ±1. The resulting readout fidelities are asymmetric: Fq = 0.890(4) for rua = 0 (limited by the 
detection efficiency and number of cycles before a spin flip) and Fi — 0.988(2) for nia = ±1 (limited by background counts 
and unwanted excitations). Because uncontrolled spin flips in the excited state decohere nearby nuclear spins, we minimize 
the number of unnecessary optical excitations by using a weak readout pulse with a maximum duration of 114 fis (~ 100 
excitations) and by switching off the laser within 2 /is (~ 2 excitations) once a photon is detected [3]. The resulting 
measurement is non-destructive: the probability that the spin prepared in nia = 0 is still in that state after a measurement 
with outcome nia = 0 is 0.992. In contrast, without dynamically stopping the laser the spin would be pumped almost 
completely to nia = ±1- For the final readout at the end of the experiment, which is allowed to be destructive, we use a 
stronger readout pulse of maximum duration 35 fis. 
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PmJ 
Pm.2 
Pm,3 



10 ), 

10 ), 

10)3 


Supplementary Figure 3: Qubit initialization. All qubits naturally start in fully mixed states pm- The ancilla is 
initialized in |0)a and a reduced SWAP operation between the ancilla and the qubit is performed, deterministically initializing 
the qubit in |0). The ancilla is then reinitialized by a 300 ps laser pulse (Reset) and the process is repeated to initialize the 
other qubits. 
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Supplementary Figure 4: Tomography sequences for three-qubit states. Examples of three-qubit expectation 
values that are measured by mapping the required correlation on the ancilla before reading it out. a, {Xi, I 2 , 13 }, b, 

{Xi, X 2 , 13 }, c, {—Xi, Y 2 , Z 3 }. Note that the phase of the last 7r/2-pulse on the ancilla depends on the number of qubits read 
out (i.e. the number of operators that are not I). The examples given here for the measurement of one- two- and three-qubit 
expectation values can be translated to any of the 63 measurements in the full three-qubit state tomography. 
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I - Prepare II - Initialize III - Encode IV - Detect V - Correct VI - Readout 


Supplementary Figure 5: Experimental sequence and logic for the QEC experiments. Example for a single round 
of quantum error correction by stabilizer measurements (as in Fig. 3). The order of the sequence is controlled in real-time by 
an ADwin microprocessor. I - The NV center is prepared in its negative charge state and on resonance with the readout and 
reset lasers (Fig. by turning on both lasers, counting the fluorescence photons and requiring a threshold to be passed (“CR 
Check”). The nuclear spin is initialized by measuring it and continuing only for outcome mi = — 1 (“MBI N”). II - The 
qubits are sequentially deterministically initialized following Fig. i III - The encoding of the logic state is a probabilistic 
process, as shown in Fig. 2a. When the wrong outcome is obtained the preparation of the experiment starts over. IV - Errors 
are detected by two stabilizer measurements. Depending on the outcome (—1 or -|-1) of each measurement, the next sequence 
to execute is communicated to the waveform generator in real time. V - Depending on which of the 4 outcomes is obtained, a 
set of gates is performed to correct errors and to map the desired expectation value onto the ancilla (Fig. |4|. VI - Finally the 
ancilla is read out. Each outcome is taken into account without post processing or post selection. Note that for the 
experiment with three rounds of error correction (two rounds of stabilizer measurements QEC, Fig. 4b), the sequence 
branches in 16 paths instead. 
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Post-selected 

results 



Deterministic 

entanglement 
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Supplementary Figure 6: Deterministic entanglement by stabilizer measurements including post-selected 
results, a, For qubits 1 and 2. b, For qubits 2 and 3. First the qubits are initialized following Fig. in |00) with fidelity 
0.878(6) for (a) and 0.910(6) for (b) (left column). Then a XX measurement is performed (Fig. Id). Depending on the 
measurement outcome feedback is applied, so that independent of the outcome the same two-qubit state is obtained, as can 
be seen by post-selecting on the two outcomes (middle column). The full result is a deterministically entangled state (right 
column). The fidelity with the desired two-qubit entangled state is 0.776(7) in (a) and 0.824(7) in (b). 
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Supplementary Figure 7: Three-qubit state tomography for |000) and the six logical states. The three qubits are 
initialized as in Fig. and encoded following Fig. 2a. The fidelities F shown are the three-qubit state fidelities and the 
shaded bars indicate the ideal states. Ideally, the 6 encoded states are all eigenstates of the XX stabilizers with eigenvalues 
-1-1, in agreement with the high values for {X\, X 2 , 13 ), {Xi, I 2 , X 3 } and {li, X 2 , X 3 } for all states and an average fidelity with 
this code subspace of 0.839(3). The logical qubit is encoded as a|0)L -|- /3|1)l, with |0)l = ^(|-|-X,-|-X,-|-X) -|- |-X,-X,-X)) and 
|1)^ = ^(|-|-X,-|-X,-|-X) — |-X,-X,-X)). The logical state expectation values are given by: {X)l = {Xi, I2,13), 

(Y)L = (Fi, Z 2 , Z 3 ), {Z)L = {Z\, Z 2 , Z 3 ) or cyclic permutations. The logical qubit fidelities for the states are 
\+X)l-. 0.916(6), \-X)l. 0.911(6), \+Y)l. 0.822(7), \-Y)l. 0.828(7), |0)i : 0.813(9) and |l)i : 0.808(9). 
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Supplementary Figure 8: Process fidelity and error syndrome probabilities for different ancilla assignments. 

Left: comparison of the process fidelities for the quantum error correction experiment in Fig. 3 for all four possible 
assignments of the ancilla states to the +1,+1 outcome of the stabilizer measurements. In Fig. 3b we average over these four 
curves. In Fig. 3c the optimal result is used (assignment {|l)a, |l)ii}). Solid lines are fits to Eq. [^taking into account Eq. 
and yield: w = 0.8(1), w = 0.71(7), w = 0.95(7) and w = 0.84(9) for the four assignments. Right: the probabilities for the 
error syndromes for each of the four ancilla state assignments. Solid lines are expected curves similar to Eqs. |17|20| based on 
the estimated initial errors in the encoded states: = 0.091(2), = 0.064(2), = 0.077(2) obtained from Fig. The 

theoretical probabilities are in good agreement with the experimental values (no free parameters). The probabilities are the 
normalized occurrences in 84000 samples for the assignments {|0)a, |0)a} and {|0)a, |l)ii} and in 28000 samples for the 
assignments {|l)a,|l)a} and {|l)a,|0)a}. 



Supplementary Figure 9: Error syndrome probabilities for naturally occurring errors. Corresponding to Fig. 4d. 
Solid lines are theoretical predictions from the unique coherence times T 2 of the individual qubits and the initial error (pin) 
determined from this data. As the stabilizer measureme nts are performed halfway the waiting time, the error probability for 
each qubit is: Pe(|) = |(1 — Exp[—(^ t^)^]). Using Eqs. 


17120 


and the measured error outcome probabilities at the hrst 

datapoint {t = 2.99 ms), we estimate the input errors at t = 0 to be = 0.049(2), p[^ = 0.0804(4) and pj® = 0.110(2). 
Dashed lines show the expected probabilities for complete dephasing. The probabilities are based on the normalized 
occurrences in 12000 samples. 
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Supplementary Figure 10: Zoom of data of Fig. 4b (a) and numerical simulation of Fig 4d (b). a, A zoom-in of 
the area of the data in Fig. 4b in which additional rounds of error-correction are advantageous, b, Numerical Monte-Carlo 
simulations for the error correction experiment of Fig 4d. The initial state is \X)l- Each qubit then coherently evolves with a 
constant detuning randomly drawn from a Gaussian probability distribution with a — \f2lT^. We take into account the 

measured longitudinal relaxation of each qubit, which approximately decays with for ancilla state |0)a (See Tab. [^, 

the error-dependent readout fidelity (Eqs.[3]|^, and a longitudinal relaxation of the ancilla with time-constant 300 ms. 

Finally, detected errors are corrected (this last step is omitted for “No feedback”). The simulation results qualitatively match 
all the main features of the observed dephasing curves (Fig. 4d). For short times, the stabilizer measurements suppress errors 
by stopping small errors from building up coherently and error correction further reduces the remaining errors. For long 
times, the stabilizer measurements halfway the sequence preferentially suppress coherent evolutions that would result in an 
error at the end of the sequence. As a result the fidelity at long times exceeds 0.5 and decays only slowly. Moreover, applying 
error correction becomes detrimental: at the moment the stabilizer measurements are applied the state is essentially random 
and no useful information about errors is extracted so that applied corrections further dephase the final state. A complete 
quantitative comparison would require detailed modelling of the full evolution of the 4-qubit system during the gates, 
initialization, stabilizer measurements, and readout sequences as well as of the longitudinal decay at short times. 
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Error probability p 




Supplementary Figure 11: Data of Fig. 3b (a) 3c (b), 4b (c) and 4d (d), without correction for the final 
readout gates, a, b, The fitted value for w is indentical as in the main text (curve shapes are not influenced by the readout 
calibration), c, Corresponding fit values for the unencoded qubit: w = —0.02(2), for 1 round: w = 0.56(6), for 2 rounds 
w = 0.64(4) and for three rounds w = 0.70(2). d, The logical qubit follows Eq. [^with T = 13.7(1) ms and n = 2.35(8). 
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